Quantum magnetism in the stripe phase: bond- versus site order. 
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It is argued that the spin dynamics in the charge-ordered stripe phase might be revealing with 
regards to the nature of the anomalous spin dynamics in cuprate superconductors. Specifically, if 
the stripes are bond ordered much of the spin fluctuation will originate in the spin sector itself, 
while site ordered stripes require the charge sector as the driving force for the strong quantum spin 
fluctuations. 
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For quite some time it has been suspected [jyg] that the 
anomalous spin dynamics of superconducting cuprates 
has to do with the 0(3) quantum non-linear sigma model 
(QNLS), describing the collective dynamics of a quantum 
anti-ferromagnet The discovery of the stripe phase 
jl opens a new perspective on these matters. Below the 
stripe-charge ordering temperature, charge fluctuations 
have to become inconsequential and the remaining spin 
dynamics should fall automatically in QNLS universal- 
ity. As will be explained, the available data suggest that 
this spin dynamics is characterized by a close proximity 
to the QNLS zero temperature transition. This enhance- 
ment of the quantum-spin fluctuations as compared to 
the half-filled antiferromagnet can have a variety of mi- 
croscopic sources. Here we will focus on the possibility 
that these are due entirely to the charge-ordering induced 
spatial anisotropy in the spin system. Although the in- 
fluence of spatial anisotropy is well understood on the 
field-theoretic level [|]](| , the charge can be bond-ordered 
or site ordered [Q and this links the spin physics of the 
stripe phase to that of coupled spin ladders IB-pi. At 
superconducting doping concentrations, bond- and site 
order translate in coupled two-leg and three-leg spin lad- 
ders, respectively. We will present an in-depth quan- 
titative analysis of both problems, showing that spatial 
spin-anisotropy has to be largely irrelevant for site-order, 
while it might well be the primary source of quantum spin 
fluctuations in the bond-ordered case. A strategy will be 
presented to disentangle these matters by experiment. 

Let us first comment on the available information re- 
garding the stripe phase spin system. The spin order- 
ing temperature appears to be strongly surpressed as 
compared to half- filling M. A first cause can be a de- 
crease of the microscopic exchange interactions. How- 
ever, the more interesting possibility is that some micro- 
scopic disordering influence has moved the antiferromag- 
net closer to the zero-temperature order-disorder transi- 
tion (quantum critical point). The few data available at 
present seem to favor the second possibility. We specif- 
ically refer to the ESR work by Kataev et. al. [[ll| on 
Lai.gg- x -yEuyGdo.oiSr x Cu04 exploiting the Gd local 



moments to probe the spin system in the CuO planes. 
Quite remarkably, little change is seen in the spin-lattice 
relaxation rate (1/Ti) at the charge ordering tempera- 
ture, T co ~ 70K. Above T co the 1/Ti is quite similar to 
that in Lai- x Sr x CuO^ where it is known from e.g. Neu- 
tron scattering that the magnetic correlation length £ is 
already quite large at the temperatures of interest: since 
the width of the incommensurate peaks is smaller than 
their separation, the correlation length is larger than the 
stripe spacing |Q. It follows that at T ~ T co a con- 
tinuum description of the spin dynamics should be sen- 
sible. Below T co 1/Ti starts to increase exponentially 
upon lowering temperature, signalling the diverging cor- 
relation length associated with the renormalized classical 
regime. Taken together, this fits quite well the expec- 
tations for a quantum antiferromagnet which is rather 
close to its quantum critical point, with a crossover tem- 
perature from the renormalized classical- to the quantum 
critical regime T* ~ T c , 
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FIG. 1. Schematic distinction between site ordered (a) and 
bond ordered (b) stripes. 

The increase of the coupling constant go, controlling 
the long wavelength fluctuations, originates in some mi- 
croscopic phenomenon. A limiting case is that charge 
can be regarded as completely static even on the scale 
of the lattice constant, such that its effect is to cause 
a spatially anisotropic distribution of exchange interac- 
tions |q||. As indicated in Fig. 1, there are two options 
|Q: the stripes can be bondov site ordered. It is expected 
that the spin dynamics associated with the hole-rich re- 
gions is characterized by a short time scale and the mag- 
netic ordering phenomena are therefore associated with 
the magnetic domains. The spin-only model of relevance 
becomes either a spin S — 1/2 Heisenberg model de- 
scribing three-leg ladders (site ordered) or two-leg lad- 
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ders (bond ordered) with uniform exchange interactions 
(J), mutually coupled by a weaker exchange interaction 
coupling (a J, a < 1). This model is explicitely, 



perature T* to the RC regime can be deduced from the 
exponential increase of the correlation length at low T, 
using 



+aJ Yl ¥i+6„, (!) 

ix=pni ,i y 

where i = (i x ,i y ) runs over a square lattice, S x = (1,0), 
5 y = (0, 1). ni measures the width of the ladder and p 
counts the ladders. 

Since the interest is in non-universal quantities as re- 
lated to the non-trivial lattice cut-off, we studied the 
model Eq.(]l|) numerically using a highly efficient loop 
algorithm Quantum Monte-Carlo method J13| supported 
with a technique of improved estimators 111]. To keep 
track of the various finite temperature cross-overs we fo- 
cussed on the temperature dependence of the staggered 
correlation length in both the directions parallel- (£ y ) and 
perpendicular (£ x ) to the stripes. We typically insisted 
on 3 * 10 4 loop updates for equilibration and (2 — 3) * 10 5 
updates for a measurement, keeping the dimensions of 
the system in the x and y directions L xy > 6 * £,x,y, 
to avoid finite size effects Jl5[. The correlation lenght 
was determined by fitting the staggered spin-spin correla- 
tion function C(f) = (— l) r * +r y (^S i+r ■ Si\ using a sym- 
metrized 2 dimensional Ornstein-Zernike form C(r) = 
A( r -V2 e -r/« + (£_ r )-i/2 e -(£-r)/£) sep arately for the x- 

(r = (r, 0),L = L x ) and y (r= (0, r), L — L y ) directions 
omitting the first few points to ensure asymptoticity. We 
checked our results against the known results for both iso- 
lated ladders by Greven et. al. [jl6| (a = 0, ni = 1,2,3) 
and the low temperature results for the isotropic (a = 1) 
limit @0- 

Since 0(3) universality is bound to apply at scales 
much larger than any lattice related cross-over scale, uni- 
versal forms for the temperature dependence of the cor- 
relation length can be used to further characterize the 
long wavelength dynamics. The absolute lattice cut-off 
is reached at a temperature (T max ) where the correla- 
tion length parallel to the stripes (£ y ) becomes of order 
of the lattice constant. However, the problem is charac- 
terized by a second cut-off: when the correlation length 
is less than the lattice constant in the direction perpen- 
dicular to the stripes (a x ), the dynamics is that of N x 
independently fluctuating spin ladders. We define T as 
the temperature where £ x ~ a x being the cross-over tem- 
perature below which the system approaches 2+1D 0(3) 
universality. In this latter regime, further cross-overs are 
present. When the effective coupling constant (go) is less 
than the critical coupling constant (g c ) a cross-over oc- 
curs from a 'high' temperature quantum critical-(QC) to 
a low temperature renormalized classical (RC) regime. 
In the QC regime £ ~ 1/T while the cross-over tem- 



e T'/T 

Z(T) « 2T*+T ' ^ 

where T* = 2irp s in terms of the spin stiffness p s (ot). 
When go > g c , the ground state is quantum disordered 
(QD) as signalled by £ becoming temperature indepen- 
dent, and the crossover temperature T' between the QC- 
and QD regimes is estimated from the approximate rela- 
tion @, 

where c y is spin wave velocity in the strong direction. 

We determined the various cross-over lines as function 
of a for the cases n; = 1,2 and 3 (anisotropic Heiscnbcrg, 
coupled two- and three leg ladders, respectively). To de- 
termine To, we used for a close to 1 the same criterium 
as for the T max determination in the isotropic problem 
(£x(T°) = 0.7 — 0.8). This becomes inconsistent for small 
a where one better incorporates the width of the ladder 
(£x(To) = ni x (0.7 — 0.8)) and we used a linear interpo- 
lation to connect smoothly both limits. We checked that 
below the To, determined in this way, both £ x and £ y 
exhibited the same dependence on temperature after an 
overall change of scale, demonstrating that the collective 
dynamics is indeed in a 2+ ID regime. 
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FIG. 2. Crossover temperatures as a function of anisotropy 
a for the coupled three leg (a) and two leg (b) spin-ladder 
models. The lines and points refer to the analytical- and nu- 
merical results, respectively, for the various scales. Notice 
that the 1-leg 'cut-off' (1D-2D cross-over) follows closely the 
results for T* . 
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In Fig. 2 we summarize our results in the form of 
a cross-over diagram as function of a and temperature, 
both for the 1- and 3-leg (Fig. 2a) and the 2-leg (Fig. 
2b) cases. Consistent with analytic predictions 0], the 
behavior is radically different for the half-integer spin 1- 
and 3 leg cases on the one hand, and the 'integer spin' 
2-leg case on the other hand. Let us first discuss the for- 
mer. Here the ground state remains in the renormalizcd 
classical regime for any finite a. The reason is obvious. 
In isolated ladders (a = 0) with an uneven number of 
legs the ground state is a Luttinger liquid exhibiting al- 
gebraic long range order and any finite ladder-to-ladder 
interaction will suffice to stabilize true long range or- 
der at T = §,0. This in turn implies a finite T* 
where the classical nature of the ground state becomes 
visible. Interestingly, our calculations indicate that T* 
and T° basically coincide for any a: at the moment the 
system discovers that it is 2+1 dimensional, the classical 
behavior sets in. Our finding that T° increases linearly 
with a for small a (Fig. 2a) confirms the scaling the- 
ory by Affleck and Halperin for this problem |l(J. The 
behavior of the spin-spin correlator for an isolated chain 
(S(x)S(O)) ~ (l/a:)exp(-£/£i), & - 1/T signals the 
approach to the Gaussian fixed point: within the ther- 
mal length £i the system exhibits algebraic long range 
order. For finite a the crossover temperature T° can be 
found using the standard mean-field consideration: at 
T° temperature becomes of order of the exchange inter- 
action between two patches of correlated spin on neigh- 
boring chains of size £i: ksT ~ a$f/^i(T°), where 
$1 = £i{T°)4>, 4> being the microscopic staggered magne- 
tization. Taking ip independent of a would yield the er- 
roneous result that T° ~ y/a. The subtlety is that when 
a is sufficiently small, the quantum dynamics within the 
correlation volume £i is already in the 2+1D regime ||. 
Using the T = result by Affleck and Halperin that 
(f) ~ y/a |[o| , we recover T° ~ a, a « 1. The other fea- 
ture worthwhile mentioning is that T° and T* are iden- 
tical for the 1- and 3-leg cases for small a's. This is in 
line with the observations by Frishmut et. al. 0] that 
these spin ladders renormalize in identical Luttinger liq- 
uids when the ladder exchange interactions are isotropic. 

In the two-leg ladders case (Fig. 2b) the quantum 
order-disorder transition occurs at a finite value of a, 
a c = 0.30(2). This is in line with the qualitative ex- 
pectations (see also [^l|,|2|). Since the isolated two leg 
ladders are incompressible spin systems, the ladder-to 
ladder interaction has to overcome the single ladder en- 
ergy gap before the two dimensional lock- in can occur. 
This critical a is rather large, and in addition, the 1+1 
D — > 2+1 D crossover temperature T° shows the upward 
curvature (T ~ y^J) previously predicted from a scaling 
analysis of the anisotropic QNLS model (AQNLS) ||. 
As a ramification, T° and T* (as well as T") separate 
and a large, genuinely 2+1 D quantum critical regime 
opens up around a c . This is in marked contrast with 



the isotropic Heisenberg model where the renormalized 
classical regime sets in essentially at the lattice cut-off 

The gross a dependences of the various cross-over tem- 
peratures can be understood by considering the AQNLS 
model obtained by taking the naive continuum limit for 
the ladder problem. An average staggered field 4> is in- 
troduced for a block of 2 x n; sites. Integrating out the 
quadratic fluctuations Q, the effective action for <j> be- 
comes the AQNLS model with anisotropic spin wave ve- 
locities, 
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3(7+3a) 
2(13+2a) 



for ni = 2 
for m = 3 

= 2 
= 3 ' 



(4) 
(5) 



where Cq is the spin wave velocity in the isotropic limit. 
The coupling constant g is a independent and the same 
as for the isotropic model. According to the scaling anal- 
ysis of Ref. the renormalized spin-stiffness becomes 
in terms of the velocities c T 



u x,yi 



Ps{a) 



c a (l-^o ) 



(6) 



where 

g c (a) = A-kaJ co/cy ^1 + -(c y aTsmh[c x /c y }/c x 

+ln[c,(l + y Jl + cl/cl)/c x /{\ + v^) 2 ])) (7) 

and p s is the spin stiffness for a — 1. According 
to Ref. ||, the cross-over scales are T* = 2irp s (a), 
T° = 2Trp s c x {g a /(4irc ) + (1 - go/g c )/c y ) and V = 
const.|p s (a)|. It turns out that for the bare coupling con- 
stant g° as determined for the isotropic case (g° = 9.1), 
the order-disorder transition occurs at a somewhat small 
value of a = 0.08, which is not surprising given the ap- 
proximations involved (one- loop level). However, by ad- 
justing go to shift a c to its numerical value (go — 11.0), 
we find a very close agreement between the numerical- 
and analytical results for the various cross-over tempera- 
tures (Fig. 2b). As can be seen from Fig. (2a), the above 
analysis also works quite well for the three-leg ladders for 
a > 0.4. Remarkably, it seems that T* switches rather 
suddenly from the AQNLS behavior at large a to the 
linear behavior expected for the Luttinger liquid regime, 
as if the topological terms start to dominate rather sud- 
denly. 

Besides its intrinsic interest, the above does have po- 
tentially important ramifications for the understanding 
of the quantum-magnetism in cuprates: bond ordering 
of stripes would imply that already at rather moder- 
ate values of the anisotropy a, spin-ladder physics alone 
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would enhance the quantum spin fluctuations substan- 
tially. This can be further illustrated by comparing the 
temperature dependence of T£ y (T) for the isotropic spin 
system a = 1 with that of the coupled two-leg ladders in 
the vicinity of the critical a (Fig. 3). This quantity can 
be directly compared with the spin-spin relaxation rate 
1 /T2G and, with some caution, also to 1 jT\ Jll],|25| (a dy- 
namical critical exponent z = 1 is only strictly obeyed in 
the QC regime). As compared to the isotropic case, the 
exponential increase of T£ (signalling the renormalized 
classical regime) is shifted to a low temperature, while 
over most of the temperature range T£(T) is constant, 
as is found in cuprates. It is noted that the 'quantum- 
critical signature' £ ~ 1/T extents in the temperature 
range above the dimensional crossover temperature T°. 
Since this regime is non-universal this should be regarded 
as a quasi-criticality. 
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FIG. 3. £ y T versus temperature for the 2-leg system, when 
the a's are close to critical point. Results for a = 0.0 (isolated 
ladders) and 1.0 (isotropic limit) are added for comparison. 
The vertical bar indicates the 1D-2D cross-over temperature. 

This is no more than suggestive. However, it points 
at a simple strategy to clear up these matters by experi- 
ments involving the static stripe phase. It should be es- 
tablished if the stripe phase is site- or bond ordered which 
can be done by NMR. Next, the a should be determined 
from neutron measurements of the spin-wave velocities, 
Eq.(||). Using these as an input, the temperature depen- 
dence of the correlation-length, as well as the NMR relax- 
ation rates, can be calculated to a high precision starting 
from a microscopic spin-only dynamics. Comparison of 
these quantities to experiment should yield insights in 
the microscopic origin of the peculiar spin dynamics in 
doped cuprates. 
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